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Abstract 

In this letter I mainly consider a finite XXZ spin chain with periodic boundary 
conditions and an odd number of sites. This system is described by Hamiltonian 

H XX z = — J^jLiWj a j+i + ^j^i+i + ^ a j a j+i}' ^ s ^ turns out, the ground state 
energy is proportional to the number of sites E = —3N/2 for a special value of the 
asymmetry parameter A = —1/2. The trigonometric polynomial q(u), the zeroes of 
which are parameters of the ground state Bethe eigenvector, is explicitly constructed. 
This polynomial of degree n = (N — l)/2 satisfies the Baxter T-Q equation. Using the 
second independent solution of this equation that corresponds to the same eigenvalue of 
the transfer matrix, it is possible to find a derivative of the ground state energy w.r.t. 
the asymmetry parameter. This derivative is closely connected with the correlation 
function < Pj&j+i >= — 1/2 + 3/2N 2 . This correlation function is related to the 
average number of spin strings for the ground state < N str i n g >= \{N — 1/N). I 
would like to stress that all the above simple formulas are not applicable to the case 
of an even number of sites which is usually considered. 

I did not care what it was all about. 
All I wanted to know was how to live in it. 

(Ernest Hemingway) 

About 30 years ago Baxter noticed that in some exceptional cases the ground state 
energy of the XYZ spin chain, which has Hamiltonian 

N 

H xyz = - i ■'■ n j a 'i- i + JyVjVj+i + J z a ]°]+\ ) > °n+i = #i , ( 1 ) 

i=i 

has the especially simple value 

lim E/N = —(J x + J y + J z ), if J x , J y and J z satisfy J x J y + J y J z + J Z J X = 0. 

in the thermodynamic limit. He later noted B, that the inversion relation gives a very simple 
eigenvalue for the 8- vertex model transfer matrix that corresponds to (0). Using standard 
notation for the Boltzmann weights of the 8-vertex model we can reformulate the Baxter's 
remark as follows. If the weights satisfy the condition 

(a 2 + ab) (b 2 + ab) = (c 2 + ab) (d 2 + ab) , (2) 
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then the transfer matrix has (up to a sign) eigenvalue T = (a + b) . Consequently Hamilto- 
nian ([[J), even for finite chains, has an exact eigenenergy of E = —N(J X + J y + J z )- A natural 
suggestion follows. Suppose that the corresponding eigenvector is the ground state vector, 
then it is probably possible to obtain interesting information for the finite chains. However, 
there are some problems. It is evident that the Baxter's remark is valid for N=l when both 
eigenvalues of the transfer matrix are T = a + b, but for N = 2, the transfer matrix does not 
have the eigenvalue T = (a + b) 2 . For simplicity let us consider the 6- vertex model which 
is the trigonometric limit (d = 0) of the 8-vertex model. Hamiltonian ([!]) is reduced to the 
XXZ Hamiltonian with a special asymmetry parameter A = —1/2: 

N 1 

H xxz = - ^{ojcr -. , + o)o]+ x - -cr|crj +1 }, a N+1 = (3) 
j'=i 

Solving Baxter's T-Q equation for N = 2 we easily find a trigonometric polynomial Q(u) of 
degree 2, but the corresponding Bethe vector does not exist. 

There are at least three ways to fix this problem and to obtain a simple eigenvalue in 
the transfer matrix spectrum. The first two have to do with modifying system (||]). 

Firstly one can modify the boundary conditions. In 1987 Alcaraz et al. || considered an 
open XXZ spin chain with a special magnetic field at the boundaries. The authors carried 
out an intensive investigation of the system. For A = — 1/2 they found a linear dependence 
between the ground state energy and the number of sites. The Hamiltonian of this chain can 
be expressed with the help of elements of the Temperlay-Lieb algebra. For A = —1/2 this 
algebra has the trivial one- dimensional representation. In the circumstances it is possible to 
explain the simplicity of the ground state energy. At present, due to multiple investigations in 
nineties, (see, for example, the paper of Hinrichsen et al. |4j] and the paper of Martin-Delgado 
and Sierra ||) this system is considered trivial. I intend to discuss all these questions in a 
future publication. Now I limit myself to a reference to paper |J, where an exact solution 
of the Baxter's T-Q equation was found. It is of importance that this solution corresponds 
to the ground state of the system. 

Secondly one can apply an additional horizontal field. This field breaks the spin reversal 
invariance of the original model. As a result, a lot of new Bethe states emerge. For a special 
value of the field strength one finds in the spectrum of the transfer matrix the simple value 
T = — (a + b) N (for even N only). The associated spin chain Hamiltonian is described by 
Perk and Schultz 0. I am indebted to Bazhanov and Baxter who informed me about this 
possibility. It was investigated in the sixth section of paper M, which is an extended version 
of paper |J. 

The two above mentioned methods deal with the trigonometric case only. I chose a third 
way. As it turns out, it is enough to consider the usual XYZ {XXZ) - spin chain 
with periodic boundary conditions, but with an ODD number of sites N — 2n+l. 
I have checked that for N = 3, 5, 7 the transfer matrix for 8-vertex model has the largest 
eigenvalue T = (a + b) N when the weights satisfy condition (0). 

Due to technical difficulties typical for 8-vertex model, I limit myself here with the trigono- 
metric case (d = 0). The existence of the above-mentioned simple eigenvalue for the ground 
state energy was discovered in the trigonometric case by Alcaraz, Barber and Batchelor [25 1. 
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In this case the 8-vertex model reduces to the 6-vertex model and formula (0) reduces to 
c 2 = a 2 + ab + b 2 . The corresponding Hamiltonian is given by The simplicity of the 
ground state in this case allows one to find simple explicit formulas for some of 
the correlations. 

Let us consider a product of z-axis spin projection operators acting at the neighbouring 
sites. One of the simplest correlations is the average of this product over the ground state 

cj(A)=<o^ +1 >. (4) 

Due to translation invariance, this correlation does not depend on j. Similar correlations 
related to the x and y axes do not depend on j also. Their values coincide, due to z-axis 
rotation invariance. Below we use the notation 

cf(£) = <ofof +1 > = <o* j0 * j+1 >. 

The simple formulae for these correlations, valid for A = — 1/2 and for odd N, 
are the main result of the paper. 

Firstly let us formulate the starting point of our calculations. It is well known that 
the Boltzmann weights of the 6-vertex model can be conveniently parametrized by spectral 
parameter u and crossing-parameter i] as 

a = sm(w + r}/2), b = sin(w — rj/2) and c = sin 77. (5) 

The asymmetry parameter is related to the crossing parameter via A = cost]. For the 
special case u = rj/2, transfer matrix T is proportional to a shift in the chain by one site 
and its eigenvalues for all states with zero momentum are sin^ r\. It is also known that 
the Hamiltonian for XX Z spin chain is related to the logarithmic derivative of the transfer 
matrix w.r.t. the spectral parameter: 

N 

H xxz = - J2i a j a j+i + ff M+i + cos V °j } = Wcos77-2sin77(-^) . (6) 

We will also use Baxter's T-Q equation 

T(u) Q(u) = sin N (u + r}/2) Q(u - r]) + sin N (u - rj/2) Q(u + rj), (7) 

where T(u) and Q(u) are eigenvalues of transfer matrix T and of auxiliary matrix Q, corre- 
sponding to a common eigenvector. In the trigonometric case 

m 

Q(u) = Y[ sm(u — Uj). (8) 

3=1 

where Uj satisfy the Bethe equation. 

Firstly , we explicitly find for an odd value of N and for a fixed value of the crossing 
parameter 77 = 2tt/3 two solutions Q(u) and P(u) of equation (0) corresponding to the 
hypothetical eigenvalue of the transfer matrix T(u) = (a + b) . We argue that the Bethe 
vector, constructed with the help of Q(u), corresponds to the ground state. Then , using the 
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result of paper M, we formulate the relation between these two solutions and the derivative 
of the largest eigenvalue of the transfer matrix T(u) w.r.t. rj. Lastly , the knowledge of this 
derivative lets us find the correlation (f|). 

Let us momentarily fix the crossing parameter rj = 2n/3 (A = —1/2) and consider 
the conjectured value T = (a + b) N . Using parametrization (^J), one can write this as 
T(u) = sin N u. Baxter's equation (|7|) for odd N takes the cyclic form 

/(«)+/(«+-£) + /(« + _) = (), (9) 

where f(u) = sin 2n+1 u Q{u). To solve the equation we follow ||. Using cross- invariance of 
the T-Q equation and the simple structure of the transfer matrix for u = rj/2, one can show 
that Q(u) is an even function and thus f(u) is an odd trigonometric polynomial of degree 
3n + 1, satisfying periodicity f(u + n) = (— l) n+1 f(u). We can therefore write 

f(u) = a\ sin(3n + l)u + sin(3n — l)u + sin(3n — 3)u . . . 

Equation (||) is satisfied if a^ u = 0. This condition implies 

n 

f(u) = ctk sin(l — 3n + 6k)u. 

k=0 

The polynomial f{u) is divided by sin 2n+1 u by definition. This fixes the coefficients It is 
clear that the first 2n derivatives have to be zeroes for u = 0. The derivatives of even order 
are zeroes trivially while the odd derivatives give 

n 

ofc(l - 3 ™ + 6A;) 2 ^ +1 = 0, n = 0, 1, . . . , n - 1. 
This system is equivalent to the condition that the relation 

n 

J2 a k {l -3n + 6fc)P((l - 3n + 6k) 2 ) = 0, (10) 

k=0 

is valid for all polynomials P(x) of degree n — 1. Let us consider n polynomials of degree 
li — l: 

n 

P r (x)= [| (x- (l-3n + 6k) 2 ). 

k=l,k^r 

Using these polynomials in formula (|T0"D we get the relations connecting the a r with a®. It 
is possible to write the answer in terms of binomial coefficients 

/(«) =f°t( n ~ h *)(iii) sin ( i - 3n + Qk >> ( n ) 

where fo is an arbitrary constant. 
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The auxiliary function g(u) = sin 2n+1 P(u) corresponding to the second independent 
solution of T-Q equation can be found by analogy 

*k«) =9°t( n i l )( n n -i) sin(2 - 3n+ 6k)u - (i2) 

One can easily check that f(u) and g(u) satisfy ODE: 

/" - 6n cot 3u f + (1 - 9n 2 )f = 0, 

g" -6ncot3ug' + {4-9n 2 )g = 0. (13) 

It happens that occasionally these equation are more convenient in calculations than explicit 
formulae (0) and flU). 

The integer m in @ is equal to the number of reversed spins and related to the z-axis 
projection of the total spin 5^ = N/2 — m. The solution Q(u) has degree m = n = (N — l)/2 
consequently the corresponding eigenvector has S 2 = 1/2. In principle we can construct 
this vector using QISM & 0. It is also known that for the antiferromagnetic XXZ chain 



with an even N the ground state has S z = [|TT| It can be analogically conjectured 
that one of the two ground states for the case of N odd has S z = 1/2. Note that in the 
interval u G (7r/3, 2n/3), the weights of the 6- vertex model ([5]) are positive. Consequently 
the components of the ground state vector are non-negative due to the Perron-Frobenius 
theorem. 



In collaboration with A. V. Razumov JT3], we have found explicit values for these com- 
ponents up to iV = 17. They are all positive. Hence, if one of the two ground state vectors 
has S z = 1/2 then we have identified it. Further we conjecture that the solution Q(u) we 
found corresponds to the ground state for iV > 17 as well. I believe that using ODE ( |l~3"D 
for f{u) it is possible to describe the distribution of the Bethe parameters and to prove this 
conjecture. 

For even N, the corresponding solution Q(u) has degree N/2 + 1 and there is no Bethe 
vector that ensure the simple eigenvalue we have discussed. 

Now we return to the main calculations. Firstly, we only know the largest eigenvalue of 
the transfer matrix T(u) = sin^w for 77 = 2tx/3. It is remarkable that the knowledge of the 
second independent solution P(u) allows us to find the derivative of this eigenvalue w.r.t. rj 
and thus the simplest correlations. 

Now we consider Baxter's T-Q equation ([?[). It can be interpreted as a discrete version 
of a second order differential equation, so we can express its coefficients in terms of two 
independent solutions [TJ, |r5fl : 



sin^ u = P{u + ri/2)Q{u - 77/2) - P(u - ri/2)Q{u + r//2), 
T(u) = P(u + rj)Q(u - rj) - P(u - rj)Q(u + rj). 

Using these relations similarly to paper ||, we can find the T-matrix derivative w.r.t rj 
^(w)|,=*r/3 = I{p{u + tt/3)Q'(u - tt/3) - P'{u + n/3)Q(u - tt/3) + 

+P(u - tt/3)Q'(u + tt/3) - P'(u - n/3)Q(u + tt/3)}. (14) 
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The derivative of the last equation w.r.t. the spectral parameter u is 

T^ u (u)\ v=2n/3 = + vr/3)Q> - tt/3) - P"{u + tt/3)Q(u - tt/3) + 

+P(w - 7r/3)Q"(« + tt/3) - P"{u - tt/3)Q(u + tt/3)). (15) 



Let us use these derivatives. Averaging (j^) over the ground state we obtain the energy 
per site that relates the correlations to the logarithmic derivative of the largest eigenvalue 
of the transfer matrix w.r.t spectral parameter u 

E ( V )/N = -2Ci - cos V Cf = cos V -^-(Tk) . (16) 

For 7] = 2n/3 and T{u) = sin^w this last equation is reduced to 

E /N = -2Ct(A = -1/2) + \ C|(A = -1/2) = -| (17) 

It is difficult to find two unknowns via one equation. However, let us differentiate equation 
(|16|) w.r.t rj. H is a Hermitian operator and so one may ignore the r\ dependence of the 
eigenvector. Hence, differentiating equation flT5|), the matrix element of the Hamiltonian, we 
can ignore rj dependence of the correlations 



cos 77 (^j + sin r]— (^) I 



U=T]/2 



Jimbo and Miwa |Td| used this method for calculation of the correlators C[ and Cf" in the 
thermodynamic limit. 

Replacing rj = 2tt/3 we obtain the last but one formula for the correlation: 



The problem is solved in principle. Due to formulas ( [H] ) and (|T5D, we can express the 
correlations in terms of the two independent solutions Q(u) and P(u) which are given (up 
to the factor sin^w ) by the formulas (|TTJ) and flT2|) . The final answer is 

< ^ > s Cf (A = -1/2) = -I + = _1 + 

Detailed calculations will be published elsewhere. 

Taking into account that the action of the operators 0"J<7j + i depends upon the relative 
orientation of the neighbouring spins, we can easily convert the last formula into a formula 
for the average number of "strings", i.e. clusters of spins with the same orientation: 

3 1 

< N stHng >= -(N — —). 
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Using equation (IT 



we get the second correlation: 
CftA = -1/2) 



< > 



8 + 8N 2 ' 



The obtained results can be generalized in different ways. Let us discuss some possibili- 
ties. 

First of all let us note that all matrix elements of the Hamiltonian ([3D are integers or 
halfintegers. The ground state energy is a halfinteger as well. It is not surprising that the 
normalization of the eigenvector can be chosen so that all its components are integers. This 
helps us to make calculations and using Mathematica we have explicitly found eigenvectors 
for iV < 17. The obtained information is contained in publication Let us mention 

only one result, related to the correlations, which are called Probabilities of Formation of 
Ferromagnetic String |T7|. All data are in agreement with the conjectured formula: 

< ai a 2 ... > (2k - 2)! (2k - 1)! (2n + fc)! (n - k)\ 



where a. 



;i+^)/2 



a k > ~" (k - 1)! (3A; - 2)! (2n - k + 1)! (n + k — 1) 



Many remarkable connections between the wave function components are noticeable. As 
it turns out, the ratio of the largest component to the smallest one is equal to the number 
of AMS (alternating sing matrices). The wonderful history of these numbers already has 
interweaved with the 6- vertex model (see, for example, JTJ|). Probably this is just the top 
of an iceberg. 

Second one can try to consider the elliptic case. 

Thirdly we can consider the inhomogeneous 6- vertex model. 

I would like to mention that the necessity to distinguish between even and odd chains 
was remarked upon by Faddeev and Takhtajan [lJJ for the XXX chain, Baake et al. for 
XXZ [ 20 1 , Bugrij |^TJ for the Ising model and so on. In papers [27| authors derived benefit 



from consideration of XXZ spin chains with an odd number of sites. As far as more recent 
results are concerned I want to mention our paper fll5|, where it was noticed that properties 



of Bethe equations depend on the parity of spin chain length. Papers of Schnack et al. 
contain numeric investigation for the spin 1/2,1,.. .,5/2 XXX chain and demonstrate the 
peculiarities of odd N. Finally I would like to mention very recent Albertini paper |28 



where author considers spin chains with open boundaries and stresses that antiferromagnetic 
quantum spin chains can in principle have very different properties according to the parity 
of their length. 
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